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Abstract —Heterogeneous networks (HetNets) with offloading is 
considered as an effective way to meet the high data rate demand 
of future wireless service. However, the offloaded users suffer 
from strong inter-tier interference, which reduces the benefits of 
offloading and is one of the main limiting factors of the system 
performance. In this paper, we investigate the use of an inter¬ 
ference nulling (IN) beamforming scheme to improve the system 
performance by carefully managing the inter-tier interference 
to the offloaded users in downlink two-tier HetNets with multi¬ 
antenna base stations. Utilizing tools from stochastic geometry, 
we derive a tractable expression for the rate coverage probability 
of the IN scheme. Then, we optimize the design parameter, i.e., 
the degrees of freedom that can be used for IN, to maximize the 
rate coverage probability. Specifically, in the asymptotic scenario 
where the rate threshold is small, by studying the order behavior 
of the rate coverage probability, we characterize the optimal 
design parameter. For the general scenario, we show some 
properties of the optimal design parameter. Finally, by numerical 
simulations, we show the IN scheme can outperform both the 
simple offloading scheme without interference management and 
the almost blank subframes scheme in 3GPP LTE, especially in 
large antenna regime. 

I. Introduction 

The next generation of wireless networks will see a signifi¬ 
cant increase in the number of wireless users and the scope of 
high data rate applications, leading to a high data rate demand. 
The conventional cellular solution, which comprises of high 
power base stations (BSs) covering a large cellular area, will 
not be able to scale with the increasing data rate demand of 
these future networks. A promising solution is the deployment 
of low power small cell nodes overlaid with high power macro- 
BSs, so called heterogeneous networks (HetNets), which are 
capable of aggressive reuse of existing spectrum assets, and 
thus lead to the support of future high data rate applications. 

Due to the larger power at the macro-BSs, more users intend 
to associate with the macro-cell tier, which causes the problem 
of load imbalancing between the macro-cell tier and the small¬ 
cell tier. Load imbalancing is recognized as one of the main 
bottlenecks for the performance of HetNets [1,2]. As such, 
to remit this problem, it is suggested to offload users from 
the heavily loaded macro-cell tier to the lightly loaded small¬ 
cell tier. Under the offloading scheme, the nearest macro-BS 
of each offloaded user, which provides the strongest signal 
among all the BSs, now becomes the dominant interferer of 
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this offloaded user. Thus, after offloading, the offloaded users 
have degraded signal-to-interference ratio (SIR). The perfor¬ 
mance of the offloaded users is one of the limiting factors 
of the network performance. Several interference management 
techniques have been considered to improve the performance 
of the offloaded users in HetNets with offloading. One such 
technique is almost blank subframes (ABS) in 3GPP LTE. In 
ABS, (time or frequency) resource is partitioned, whereby the 
offloaded users and the other users are served using different 
portions of the resource. The performance of ABS in HetNets 
with offloading was analyzed in [3] using tools from stochastic 
geometry. Another interference management technique was 
proposed in [4] to reduce the interference to offloaded users, 
where each offloaded user is jointly served by both of its 
nearest macro-BS and nearest pico-BS. The performance of 
this joint transmission scheme in large HetNets was analyzed 
using similar tools as in [3]. Note that these works only 
considered HetNets with single-antenna BSs. 

Deploying multiple antennas in HetNets is another solution 
to provide high data rates for future wireless service. Futher- 
more, with multiple antennas, more effective interference man¬ 
agement techniques (e.g., coordinated beamforming [5]) can 
be implemented. For example, a HetNet with a single multi¬ 
antenna macro-BS was investigated in [6,7], where the multi¬ 
ple antennas at the macro-BS are used for serving its scheduled 
users as well as mitigating interference to the receivers in small 
cells using different interference coordination schemes. These 
schemes have been analyzed and shown to have performance 
improvement. However, since only one macro-BS is consid¬ 
ered, the results obtained in [6,7] can not reflect the macro¬ 
tier interference, and thus can not offer accurate insights for 
practical networks. In [8], interference coordination among a 
fixed number of neighboring BSs was investigated in downlink 
large multi-antenna HetNets. However, this scheme may not 
fully exploit the spatial properties of the interference in large 
HetNets. Moreover, offloading was not considered in [8]. So 
far, it is still not clear how the interference coordination 
schemes and the system parameters affect the performance 
of large multi-antenna HetNets with offloading. 

In this paper, we consider offloading in downlink two- 
tier large stochastic HetNets with multi-antenna BSs, and 
investigate the use of interference nulling (IN) beamforming to 
improve the performance of the offloaded users. The scheme 
has a design parameter, which is the degree of freedom U 


that can be used at each macro-BS for IN to its offloaded 
users. In particular, each macro-BS utilizes the low-complexity 
zero-forcing beamforming (ZFBF) precoder to suppress inter¬ 
ference to at most U offloaded users as well as boost the 
desired signal to its scheduled user. Note that interference 
coordination using beamforming technique in large stochastic 
HetNets causes spatial dependency among macro and pico 
cells, which is quite difficult to analyze in general [7]. In this 
paper, by adopting appropriate approximations and utilizing 
tools from stochastic geometry, we first present a tractable 
expression for the rate coverage probability of the IN scheme. 
To our best knowledge, this is the first work analyzing the 
interference coordination technique in large stochastic multi¬ 
antenna HetNets with offloading. To further improve the rate 
coverage probability, we consider the optimization of the 
design parameter. Note that optimization problems in large 
HetNets with single-antenna BSs were investigated (see e.g., 
[9]). In [9], the objective function is relatively simple, and 
bounds of the objective function are utilized to obtain near- 
optimal solutions. The optimization problem in large multi¬ 
antenna HetNets we consider is an integer programming 
problem with a very complicated objective function. Hence, 
it is quite challenging to obtain the optimal solution. First, for 
the asymptotic scenario where the rate threshold is small, by 
studying the order behavior of the rate coverage probability, we 
prove that the optimal design parameter converges to a fixed 
value, which is only related to the antenna number difference 
between each macro-BS and each pico-BS. Next, for the gen¬ 
eral scenario, we show that the optimal design parameter also 
depends on other system parameters. Finally, by numerical 
simulations, we show the IN scheme can outperform both the 
simple offloading scheme without interference management 
and ABS in 3GPP LTE, especially in large antenna regime. 

II. System Model 

A. Downlink Two-Tier Heterogeneous Networks 

We consider a two-tier HetNet where a macro-cell tier is 
overlaid with a pico-cell tier, as shown in Fig. 1(a). The 
locations of the macro-cell BSs and the pico-cell BSs are 
spatially distributed as two independent Homogeneous Poisson 
point processes (PPPs) Tq and <f >2 with densities Ai and A 2 , 
respectively. The locations of the users are also distributed 
as an independent homogeneous PPP <F n with density A u . 
Without loss of generality (w.l.o.g.), denote the macro-cell 
tier as the 1st tier and the pico-cell tier as the 2nd tier. We 
consider the downlink transmission. Each macro-BS has N\ 
antennas with a total transmission power Pi, each pico-BS 
has N 2 antennas with a total transmission power P 2 (< Pi), 
and each user is equipped with a single antenna. We consider 
both large-scale fading and small-scale fading. Specifically, 
due to large-scale fading, in the jth tier (j = 1, 2), transmitted 
signals with distance r are attenuated by a factor ^4j, where 
OLj > 2 is the path loss exponent of the jth tier. For small- 
scale fading, we assume Rayleigh fading channels, i.e., each 
element of channel vectors is distributed as CA/"(0,1). 


B. User Association 

We assume open access [1]. Due to the larger power 
at the macro-BSs, the load imbalancing problem arises if 
the user association is according to the long-term average 
received power (RP). To remit the load imbalancing problem, 
the bias factor Bj (j = 1,2) is introduced to each tier, 
where P 2 > Pi, to offload users from the macro-cell tier 
to the pico-cell tier. Specifically, user i (denoted as uf) is 
associated with the BS which provides the maximum long¬ 
term averaged biased-received-power (BRP) (among all the 
macro-BSs and pico-BSs) at the user. Here, the long-term 
averaged BRP is defined as the average RP multiplied by a 
bias factor. This associated BS is called the serving BS of 
user i. Note that within each tier, the nearest BS to user i 
provides the strongest long-term averaged BRP in this tier. 
User i is thus associated with the nearest BS in the j*th 
tier if 1 j* = arg max je ^ 12 }PjBjZ i ™\ where Z^j is the 
distance between user i and its nearest BS in the j th tier. 
From the criterion described above, we observe that, for given 
{Pj}’ {^ij} and {ay}, the user association is only affected 
by the bias factor ratio. Thus, w.l.o.g., we assume B\ = 1 and 
P 2 = B > 1. After user association, each BS schedules its 
associated users according to TDMA, i.e., scheduling one user 
in each time slot, so that there is no intra-cell interference. 

According to the above mentioned user association pol¬ 
icy and the offloading strategy, as illustrated in Fig. 
1(b), all the users can be partitioned into the follow¬ 
ing three disjoint user sets: i) set of macro-users : U\ = 
{ Ui\P\Zff 1 > BP 2 Zf 2 2 }, ii) set of unoffloadedpico-users : 
U 2 o = { u i\P2Zf 2 2 > Pi^i* 1 }’ iii) set of offloaded users: 
U 20 = {ui\P 2 Z -? 2 <PiZ~p <BP 2 Z~? 2 }, where the 
macro-users are associated with the maco-BSs, the unoffloaded 
pico-users are associated with the pico-BSs (even without 
bias), and the offloaded users are offloaded from the macro- 
BSs to the pico-BSs (due to bias B > 1). Moreover, U2 = 
U 2 o U ^20 is the set of pico-users. 

C. Performance Metric 

In this paper, we study the performance of the typical user 
denoted as 2 uq , which is located at the origin and is scheduled. 
Since HetNets are interference-limited, in this paper, we ignore 
the thermal noise in the analysis. Note that the analysis 
including thermal noise can be obtained in a similar way. Let 
Ro = log 2 (1 + SIRo) denote the rate of the typical user, 
where W is the available resource (e.g., time or frequency), 
Lo is the total number of associated users (i.e., load ) of the 
typical user’s serving BS, and SIRo is the SIR of the typical 
user. We investigate the rate coverage probability of the typical 
user, which is defined as the probability that the rate of the 
typical user is larger than a threshold [3], i.e., 

K(t) = Pr (Ro > r) = Pr (A log 2 (1 + SIRo) > r) (1) 

1 In the cell selection procedure, the first antenna is normally used to 
transmit signal (using the total transmission power of each BS) for BRP 
determination according to LTE standards. 

2 The index of the typical user and its serving BS is denoted as 0. 
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(b) User Set Illustration 

Fig. 1 . System model (U = 1 ) and user set illustration. 


where r is the rate threshold. The rate coverage probability is 
a desired performance metric for applications with strict rate 
requirements, e.g., video services [ 2 ]. 


III. INTER-TIER INTERFERENCE NULLING 

In HetNets with offloading, the offloaded users (to the 
pico-cell tier) normally suffer from stronger interference than 
the macro-users and unoffloaded pico-users. 3 The dominant 
interference to each offloaded user, which is caused by its 
nearest macro-BS, is one of the limiting factors of the system 
performance. In this section, we first investigate an inter-tier 
IN scheme to improve the performance of these offloaded 
users. Then, we obtain the SIR under the IN scheme. 


A. Strategy Description 

We consider an inter-tier IN scheme to mitigate the dom¬ 
inant interference to the offloaded users. We first partition 
the offloaded users U20 into two sets. Note that all the 
offloaded users may not be scheduled by their nearest pico- 
BSs simultaneously, as each BS schedules one user in each 
time slot. We refer to the offloaded users, which are offloaded 
from a macro-BS and are scheduled by their nearest pico- 
BSs in a particular time slot, as active offloaded users of 
this macro-BS. In this scheme, each macro-BS conducts IN 
to some of its active offloaded users in a particular time slot, 
which are referred to as IN offloaded users of this macro-BS. 
We refer to the remaining offloaded users of each macro-BS 
as non-IN offloaded users of this macro-BS. Under the IN 
scheme, in a particular time slot, as illustrated in Fig. 1 (b), 
the offloaded users U20 can be divided into two sets, i.e., 
U20 = U20C U^2oc> where U20C is the set of IN offloaded 
users and U 2 oc i s th e set °f non-IN offloaded users. 

We now discuss how to determine the IN offloaded users 
of each macro-BS. In particular, each macro-BS, which has 

3 For each offloaded user, its nearest macro-BS, which provides the strongest 
long-term averaged RP, now becomes the dominant interferer of this offloaded 
user. However, for each macro-user or unoffloaded pico-user, the BS which 
provides the strongest long-term averaged RP is its serving BS. Therefore, 
the offloaded users suffer from the strongest interference. 
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Ni transmit antennas, makes use of at most U (U < N±) 
DoF to conduct IN. Note that U is the design parameter of 
this scheme. When U = 0, the IN scheme reduces to the 
one without interference management technique. Specifically, 
let U 20 a ,£ denote the number of active offloaded users of 
macro-BS £, each of which is scheduled by a different pico- 
BS. If U20 a ,t < U, macro-BS £ performs IN to all of its 
U 20 a ,i active offloaded users using U2o a ,£ DoF. However, if 
U 20 a ,i > U, macro-BS I randomly selects U out of U 2 o a ,i 
active offloaded users according to the uniform distribution to 
perform IN using U DoF. Hence, macro-BS £ performs IN 
to U 2 oc,i = min ([/, U20 a /) out °f U 2 o a ,e active offloaded 
users. Note that the DoF used for IN (called IN DoF) at macro- 
BS £ is U 2 oc,i- The remaining N\ —U 2 oc,e DoF at macro-BS 
£ is used for boosting the signal to its scheduled user. 

Now, we introduce the precoding vectors at macro-BSs 
and pico-BSs, respectively. First, each macro-BS utilizes the 
low-complexity ZFBF vector to serve its scheduled user and 
simultaneously perform IN to its IN offloaded users. Specif- 
ically, denote = [h^ g Mi g M2 • • • gi,te 200 ,<] > 

where 4 h 1:i ~ CAT N i,i(Oat iX i, IjvJ is the channel vec¬ 
tor between macro-BS £ and its scheduled user, and 
gi ~ CJ\Tn 1% i (O^xi? I/Vi) denotes the channel vector 
between macro-BS £ and its IN offloaded user i (i = 
1,..., U 2 oc,i)- The precoding matrix is designed to be 

Wi^ = . Then, the ZFBF vector at 

macro-BS £ is designed to be fi^ = ^ 7 ^ where wi^ 
is the first column of the precoding matrix Wi^. Next, 
the beamforming vector at each pico-BS is utilized to serve 
its scheduled user. Specifically, we use the maximum-ratio 
transmission scheme. The beamforming vector at pico-BS £ 
is h.t = where h 2 ,<? ~ CNn 2 ,i (0jv 2 xi, Tv 2 ) is the 

channel vector between pico-BS £ and its scheduled user. 

B. SIR of the Typical User 

We now obtain the SIR expressions of the typical user 
Wo e U k , for all k G JC = {1,20,20(7,20(7}. Specif- 
ically, the SIR of the typical user uo G Uk is given by 5 

p 0 k | h rt f. I 2 

I*b'fc ’ 0 3 k’°\ 

SIRfc.o = -—--p - ’ where is the 

^ =1 W( Ai)XBfc | D .^ o p- 

channel vector between uq G Uk and its serving BS, f Jfc? o is 
the beamforming vector at the serving BS of ^0 G Uk, with 

KoU,o| ~ Gamma (M/c, 1 ), jk and Mk given by Table 
I, is the channel vector between BS £ in the jth tier 

4 The notation X A Y means that X is distributed as Y. 

5 In this paper, all maco-BSs and pico-BSs are assumed to be active. The 
same assumption can also be seen in existing papers (see e.g., [2, 3 ]). 
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uo G U k and its serving BS, \Dj^\ is the distance between 
BS £ in the jth tier and uo, and B k (given by Table I) is the 
BS set which does not cause interference to uo G li k . Here, 
Bj ,o is the nearest BS of uo in the jth tier. 

IV. Rate Coverage Probability Analysis 

In this section, we investigate the rate coverage probability 
of the IN scheme. We notice that the rate coverage probability 
of the IN scheme is dependent on the probability mass 
functions (p.m.f.) of the active offloaded users U 2 o a ,o °f the 
serving macro-BS of uq when uo G U\ and the p.m.f. of the 
active offloaded users U 2 O a ,0 of the nearest macro-BS of uo 
when uq G U 20 • Note that the p.m.f.s of U 2 o a ,o and U 2 O a ,0 
are related to Poisson distribution. But the exact distributions 
are unknown, as it is difficult to calculate due to the coupling 
between the active offloaded users and the pico-BSs. Here, 
we address this challenge by first deriving approximations for 
p.m.f.s of U 2 o a ,0 and U 2 O a , 0 • Based on the approximations 
and according to total probability theorem, we have the rate 
coverage probability 7 Z(r) as follows: 6 

Theorem 1: The rate coverage probability is 



Fig. 2. Rate coverage probability vs. rate threshold r for different bias 
factors B, and with a\ = 0.2 = 4, ^ = 10 dB, Ah = 8, N 2 = 4, U = 4, 
W = 10 MHz, Ai = 0.0001 nodes/m 2 , and A 2 = 0.0005 nodes/m 2 . 
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Note that the expression in (2) is difficult to compute and 
analyze due to the infinite summations over n. To simplify 
the rate coverage probability expression 7 Z(r) in (2), we use 
the mean of the random load (i.e., E [To,i] or E [£ 0 , 2 ]) to 
approximate it (i.e., Lo,i or £ 0 , 2 ) [3]. The simplification is 
thus achieved due to the elimination of the infinite summation 
over n. The rate coverage probability with mean load approx¬ 
imation (MLA) is given by: 

Corollary 1: The rate coverage probability with MLA is 

H(t) =AiHi(t) + A 2 d'R- 2 d ( T ) + +oPr (£ 200 , 0 ) ^ 2 oc(r) 
+ A 2 o (1 - Pr (£ 2 oc,o)) ’k 2 oc( T ) ( 12 ) 

where Hi(t) = + (/ ( E[Z + lT )), ^ 2 o( T ) = 
S 2 o (/(t 1 )), ^2 oc(t) = 5 20 c(/(%^)), 

and H 2 qc(t) = S 20C (/ ( E[L + ]r )) with S k (-) given 
by (7)—(10). Here, E [L 0 ,i] = 1 + 1.28 ^a [3], and 

E[Z 0>2 ] =1 + 1.28At 2 - ' 

Fig. 2 plots the rate coverage probability vs. rate threshold 
r for different B. We see that the Analytical’ curves (i.e., 
7 Z(t) in (2)) closely match with the ‘Monte Carlo’ curves, 
although 7 Z(r) is derived based on the p.m.f. approximations 
of U 2 O a ,0 and U 2 O a , 0 • Moreover, we observe that the rate 
coverage probability with MLA (i.e., 1Z(r) in (12)) provides 
sufficient accuracy, especially when r is not very large. Hence, 
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for analytical tractability, we will investigate the rate coverage 
probability with MLA H(t) in the remaining part of this paper. 

V. Rate Coverage Probability Optimization 

In this section, we present some properties of the optimal 
design parameter which maximizes the rate coverage probabil¬ 
ity for small rate threshold regime and general rate threshold 
regime, respectively. 

For a fixed B , the optimal design parameter 17* (r) is 
defined as follows: 7 

U*(t) = arg max 72,(17, r) . (13) 

U E{0,l,...,IVi —1} 

Note that (13) is an integer programming problem. Moreover, 
the objective function 72,(17, r) in (12) is very complicated. 
It is thus quite challenging to obtain the closed-form optimal 
solution to the problem in (13). 

Now, we address this challenge and characterize the opti¬ 
mality property of 17* (r). Let A72,(17, r) = 72,(17, r) — 72,(17 — 
1 , r) denote the rate coverage probability change when the 
design parameter is changed from U — 1 to U. We can show 

AK{U,t) = A 2 oA-JZ 2 o(U,t) -Ai |A7£i(i7, t)| (14) 

where A£i(17,t) = Ki(U,t) - H^U - l, t) < 0 de- 
notes the rate coverage probability change of a macro¬ 
user, and AH 20 {U,t) = 7 Z 2 o{U,r) - 72, 2 o{U -l ,r) > 
0 denotes the rate coverage probability change of an of¬ 
floaded user. Here 8 , 1Z2o{U,r) = Pr (17, £ 200 , 0 ) P 20 c(r) + 
(1 — Pr (17, £ 200 , 0 )) ^ 20 c( r )- We observe from (14) that 
when U is increased by 1, A72,(17, r) can be decomposed 
into two parts, namely, the “gain” A 2 oA1Z 2 o(U,t) and the 
“penalty” A\ |A72,i(17, r)|. Whether A 72,(17, r) is positive or 
not depends on whether the “gain” dominates the “penalty” 
or not. In the following, to maximize 72,(17, r), we study 


the properties of A 72,(17, r) in (14) w.r.t. U by comparing 
A2oA112o(U,t) and Ai |A72,i(17, r)|. 

1) Rate coverage probability optimization when r 0: 
We first characterize the asymptotic behavior of A 1Z 2 o(U,t) 
and |A72i(?7, t)| when r -A 0, which is shown as follows: 

Lemma 1: When r -A 0, we have 9 AH 2 o(U ,t) = 
0 (r^ 2 ) and | A72i(?7, r) | = 0 (t n ^~ u ). 

According to (14), Lemma 7, and noting that A20 and A\ 
are independent of r, we have 


6 ( 7 ^) > 0 , u <n 1 -n 2 

AK(U,t) = 0 (r^ 2 ) - 0 (r^ 2 ) , U = TVi - N 2 • (15) 

0 (r Nl ~ u ) <0, U >N x -N 2 


Since 17* (r) satisfies A72,(17*(r), r) > 0 and A 72,(17* (r) + 
l,r) < 0, we see from (15) that U*(r) should be in the set 
{Ni — N 2 — 1,Ni — N 2 }, and the exact value of U* (r) depends 
on whether A 72,(17, t) is positive or not when U = N\ — N 2 
(i.e., the second case in (15)), i.e., whether the coefficient in 
0 ( t n 2 ) of A 2 oAlZ 2 o{U , r) (i.e., the first one) is larger than 
that in 0 (■ t n2 ) of A\ \AHi(U,t)\ (i.e., the second one) or 
not. According to the above discussions, the optimal design 
parameter U*(r) is given in the following theorem: 

Theorem 2: When r 0, the optimal design parameter 
U*(t) US, Where US E {A^i -N 2 -1,N 1 - N 2 }. 

Theorem 2 shows that when r -A 0, the optimal design 
parameter U*(t) converges to a fixed value in the set {TVi — 
N 2 — l,7Vi — N 2 }, which is only related to the number of 
antennas at each macro-BS and each pico-BS. Fig. 3(a) plots 
the optimal design parameter U*(r) vs. rate threshold r for 
different B. We can see that U*(t ) converges to a fixed value 
US E {TVi — N 2 — 1, TVj — N 2 } when r is sufficiently small. 
Specifically, we can see that 17* (r) = Ni — N 2 = 3 at B = 4.6 
dB, and 17*(r) = N ± - N 2 - 1 = 2 at B = 2.5 dB. These 
observations verify Theorem 2. In addition, we see that 17* (r) 


7 We make explicit the dependence of 7 Z(r) on U in this section. 

8 Here, we make explicit the dependence of Pr (£200,0) on U. 9 f( x ) — © 0(®)) means that lim^^o = c where 0 < c < 00. 
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Fig. 3. Optimal design parameter U*(t), and with a\ = a 2 = 3, ^ = 10 
dB, W = 10 x 10 6 Hz, N! = 5, N 2 = 2, X u = 0.01 nodes/m 2 ! Ai = 
0.0001 nodes/m 2 , and A 2 = 0.0015 nodes/m 2 . 
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(a) N± = 8 , N 2 = 6 , r 7 *(t) = 0.01 (b) Ah = 18, N 2 = 16, 77 * (r) = 
at ^ABS °' 19 at ^ABS 

Fig. 4. Rate coverage probability vs. bias factors B, and with a\ = 4.5, 
cu 2 = 4.7, = 13 dB, W = 10x 10 6 Hz, r = 5xl0 5 bps, Ai = 0.00008 

nodes/m 2 , A 2 = 0.001 nodes/m 2 , and X u = 0.05 nodes/m 2 . In the figures 
on top, the points at Bj* n , Bf =Q , and Bg are highlighted using black 
ellipse. In the hgures at the bottom, the rate coverage probability of each user 
type in different schemes are plotted at Bj* n , B^_ Qt and F?XbS’ respectively. 

is larger for a larger B. 

2) Rate coverage probability for general r: Unlike the case 
for small r, the “gain” A 2 oA 1 Z 2 o(U,t) and the “penalty” 
Ai\A1Zi(U,t)\ are not determined by the order terms 
0 ( t n2 ) and 0 ( r N i ~ u ) ? respectively. The optimal design 
parameter U*(r) thus also depends on other system param¬ 
eters besides Ni and N 2 . Fig. 3(b) plots the rate coverage 
probability with ML A vs. U for different B. We can see that 
besides N\ — N 2 — 1 and N 1 —N 2 , the optimal design parameter 
U*(r) can also take other values in set {0,1,..., Ni — 1}. In 
particular, from Fig. 3(b), we see that U*(r) can be 0 (at 
B = 2 dB), 2 (at B = 5 dB), and N x - 1 = 4 (at B = 10 
dB). Interestingly, Fig. 3(b) indicates that, for general r, the 
optimal parameter U*(r) increases with B , which is consistent 
with the case for small r in Fig. 3(a). 

VI. Numerical Results 

In this section, we compare the IN scheme under U*(r) with 
the simple offloading scheme without interference management 
(i.e., U — 0) and multi-antenna ABS under 77 * (r) (i.e., the 
multi-antenna version of ABS in [3]). Specifically, 77 G (0,1) 
is the continuous control variable for ABS, where 77 represents 
the resource fraction for serving the offloaded users only, 
while 1 — 77 represents the resource fraction for serving the 
macro-users and unoffloaded pico-users simultaneously. In the 


simulation, the optimal control variable 77 * (r) of ABS is 
obtained by bisection method with Ni iterations, while the 
optimal (discrete) design parameter 17* (t) of the IN scheme 
is obtained by exhaustive search over {0,1,..., Ai — 1}. 

Fig. 4 plots the rate coverage probability vs. bias factor B 
for all the three schemes. Note that under the parameters used 
in the simulation, we have sufficient offloaded users, and the 
dominant macro-interferer is sufficiently strong compared to 
the remaining macro-interferers. We see from Fig. 4 that the 
optimal rate coverage probability (maximized over B) of the 
IN scheme is larger than those of both the simple offloading 
scheme and ABS. 10 We denote the optimal B for the IN 
scheme, simple offloading scheme, and ABS as F>j* n , Bf =0 , 
and L?ab S , respectively. In particular, when B^ N is sufficiently 
large (i.e., sufficient offloaded users), if Ni is sufficiently large 
(e.g., Ni = 8 ), the IN scheme can achieve good performance 
gains over both the simple offloading scheme and ABS, as 
the penalty of sacrificing some DoF for IN becomes minor. 
From the figures at the bottom of Fig. 4(a) and Fig. 4(b), we 
see that the rate coverage probability of the offloaded users 
in the IN scheme at B± N is larger than those in the simple 
offloading scheme at Bf =0 and in ABS at F>Xbs- This is 
because the offloaded users in the IN scheme do not have (time 
or frequency) resource sacrifice, and their dominant macro- 
interferers can be cancelled. However, the offloaded users 
in ABS suffer from resource limitations, and the offloaded 
users in the simple offloading scheme suffer from the strong 
interference caused by their dominant macro-interferers. 

References 

[1] H. S. Jo, Y. J. Sang, R Xia, and J. G. Andrews, “Heterogeneous cellular 
networks with flexible cell association: A comprehensive downlink SINR 
analysis,” IEEE Trans. Wireless Commun., vol. 11, no. 10, pp. 3484- 
3495, Oct. 2012. 

[2] S. Singh, H. S. Dhillon, and J. G. Andrews, “Offloading in heterogeneous 
networks: modeling, analysis, and design insights,” IEEE Trans. Wireless 
Commun., vol. 12, no. 5, pp. 2484-2497, Mar. 2013. 

[3] S. Singh and J. G. Andrews, “Joint resource partitioning and offloading 
in heterogeneous cellular networks,” IEEE Trans. Wireless Commun., 
vol. 13, no. 2, pp. 888-901, Feb. 2014. 

[4] A. H. Sakr and E. Hossain, “Location-aware cross-tier coordinated 
multipoint transmission in two-tier cellular networks,” to appear in IEEE 
Trans. Wireless. Commun. 

[5] Y. Cui, V. K. N. Lau, and Y. Wu, “Delay-aware BS discontinuous 
transmission control and user scheduling for energy harvesting downlink 
coordinated MIMO systems,” IEEE Trans. Signal Process., vol. 60, 
no. 7, pp. 3786-3795, Jul. 2012. 

[6] K. Hosseini, J. Hoydis, S. t. Brink, and M. Debbah, “Massive MIMO and 
small cells: How to densify heterogeneous networks,” in Proc. of IEEE 
Int. Conf. on Commun. (ICC), Budapest, Jun. 2013, pp. 5442-5447. 

[7] A. Adhikary, H. S. Dhillon, and G. Caire, “Massive-MIMO meets 
HetNet: Interference coordination through spatial blanking,” submitted 
to IEEE J. Select. Areas Commun. 

[8] P. Xia, C. H. Liu, and J. G. Andrews, “Dow nlin k coordinated multi¬ 
point with overhead modeling in heterogeneous cellular networks,” IEEE 
Trans. Wireless Commun., vol. 12, no. 8, pp. 4025-4037, Aug. 2013. 

[9] Y. Lin and W. Yu, “Joint spectrum partition and user association in multi¬ 
tier heterogeneous networks,” in Proc. of Conference on Information 
Science and Systems (CISS), Princeton, NJ, Mar. 2014, pp. 1-6. 

[10] Y. Wu, Y. Cui, and B. Clerckx, “Analysis and optimization of inter-tier 
interference coordination in downlink multi-antenna HetNets with 
offloading,” 2014. [Online]. Available: http://arxiv.org/abs/1411.3271 


10 Note that the IN scheme may not provide gains in all scenarios, as 
suggested in Fig. 3. 



























